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Abstract 

The application of the divergence theorem in a non-coordinated basis is 
shown to lead to a corrected variational principle for Class B Bianchi cos- 
mological models. This variational principle is used to construct a Hamil- 
tonian formulation for diagonal and symmetric vacuum Class B models. 
[Note: This is an unpublished paper from 1984 that might be useful to 
anyone interested in Class B Bianchi models]. 



* Based in part on a thesis by S. Waller submitted to the University of Texas at Austin in partial 
fulfilment of the requirements for the Ph. D. degree. 
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I. INTRODUCTION 



For the past decade there has been a low-level but continuing interest in Bianchi-type 
cosmological models which admit a cosmic time, a gauge in which — 0, and whose t = 
const, surfaces are homogeneous three-spaces of Bianchi types I-IX JTJ. This interest, for 
the most part, centers on the mathematical structure of the Einstein equations for these 
models instead of on astrophysical applications, since these models could only represent an 
unobservable portion of the actual universe (although a search for observable consequences 
of this primitive stage continues). The advantage of the Bianchi-type models is that they 
provide a model of general relativity which has true dynamics, even though the metric 
components in the invariant basis are fuctions of time only. These models allow one to 
study problems in general relativity in a situation where they are soluble exactly (or at 
least qualitatively). 

One line of study has been the Hamiltonian formulation of the Einstein equations for 
these models which was begun with an eye toward approximate solutions and the con- 
struction of a model quantum theory for general relativity. The Hamiltonian formulation 
for Types I and IX based on the Arnowitt-Deser Misner (ADM) formulation of general 
relativity is due to Misner ||, and was extended to all models of Ellis-MacCallum |4[] 
Class A by Ryan ||. In Ref. || an extension was made to Class B models, but MacCal- 
lum and Taub || pointed out that the naive application of the ADM formalism yields a 
Hamiltonian formulation that gives incorrect Einstein equations for these models. This 
pathology excited some interest, and a number of corrected variational principles have 
been proposed for these models |7|,|8|,|9|JT0|| 

In this paper we plan to make a detailed review of the problem of the Hamiltonian 
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formulation of Class B models from a slightly different point of view. All of the variational 



principles in Refs. [7-10] have what Jantzen JT0| calls "non-potential" terms in their 
variational principles, and there is at least one other problem that these models have that 
makes the study of this Hamiltonian formulation interesting. We plan to emphasize that 
the non-potential terms depend on the treatment of certain terms in the Einstein action, 
and that no extra correction terms are needed. In fact, with the proper treatment of 
the divergence theorem in a non-coordinated basis, the correction becomes obvious. This 



correction provides a foundation for the idea of one of us JTIJ that the problem of Class B 
models is due to problems of variational principles in a non-coordinated basis, although 
not in the form originally proposed. 

A model problem for Class B models that suffers from all the same problems (although 
in a slightly different form) is that of conformal metrics in general relativity. Consider 
the metric 

ds 2 = e 2X{x ' t \ri pLP dx»dx v ). (1.1) 

The Einstein tensor is 

= 2d^d v X - dpXdvX - 2 Vflu ( v a ?d a dpX + ^n aP d a Xd p X). (1.2) 
The Einstein action in vacuum becomes 

1 = ToV / ^~9 Rd * x = ~tJ e2X [ Q V al3 d a dpX + 4r] al3 d a XdpX}d 4 x. (1.3) 

Here it is easy to see one of the pathologies we will see in Class B models: Reduced actions 
may lack equations. Varying (1.3) with respect to A gives only one equation, and setting 
G^u from (1.2) equal to zero gives nine more, and only in cases where these nine are 
automatically zero or redundant (as in k = FRW models) will (1.3) yield the correct 
Einstein equations. 



The second difficulty that we will encounter in Class B models can be illustrated by 
rewriting the action (1.3) in an orthonormal basis, = e~ x (d/dx tl ), uj^ = e x dx fl . We 
will write the action of on a function A as e^A = A,^ = e~ x d ii A. The volume element 
d A x = e~ AX u° A u 1 A uj 2 A uj 3 , so the action becomes 

I = — [ [Qrj^X a p + 10r] a/3 X a X pP° A uj 1 A uo 2 A uo 3 . (1.4) 

07T J 

If one treats the derivatives e M like ordinary derivatives and varies (1.4), one finds as the 
vacuum equations 

t P \a,P = = V af3 d a dp\ - TJ^daXdpX. (1.5) 

Varying (1.3) one finds 

Gv^dadpX + 4r] Qf3 d a \dp\. 

The obvious problem is treating e M derivatives as though they were partial derivatives 
in integration by parts. This is the second problem encountered in Class B models: 
Integration by parts in non- coordinated bases can lead to incorrect equations. In Sec. II 
we will show that the answer to this problem lies in the proper use of the divergence 
theorem in a non-coordinated basis based on the model of Spivak []T2 or Lovelock and 
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The main thrust of this paper is to show how the application of the correct divergence 
thoerem can help in the understanding of the non-potential terms in the Hamiltonian for- 
mulation of Class B models. In fact, the only non-potential terms are those that arise from 
the derivatives of the connection coefficients in the Ricci scalar in the Einstein variational 



principle. As a practical excercise we apply the results to diagonal and symmetric [14] 
Class B metrics. While we feel that our presentation is more didactic, some of the results 
for diagonal and symmetric models can be found in the exhaustive review of Jantzen |K|. 



Note that Jantzen manages to achieve a Hamiltonian formulation without non-potential 
terms valid for certain metric variables, at the cost (admmittedly slight) of introducing 
g 0i terms which must be found by integrating supplementary equations. 

The plan of the paper is as follows. In Sec. II we discuss the divergence theorem in 
a non-coordinated basis and its application to the Hamiltonian formulation of Class B 
models. In Sec. Ill we consider vacuum Class B models with diagonal space metrics, 
showing which of them allow solutions. Certain models do not, and in Sec. IV we show 
that all models with symmetric metrics do allow vacuum solutions. 



II. THE DIVERGENCE THEOREM IN A NON-COORDINATED BASIS AND 
THE VARIATIONAL PRINCIPLE FOR CLASS B MODELS. 



As mentioned in the introduction the corrected variational principles that have been 
proposed for Class B models [7-10] have non-potential terms (with the exception of the 



final Hamiltonian of Jantzen []I0|). The first calculation of correction terms was due to 
Sneddon J7| , who called, the non-potential terms "surface terms" . Our approach is slightly 
different, and, we feel, somewhat more didactic. 

The basic cause of the problem for Class B models is that the divergence theorem in 



a non-coordinated basis [12,13] does not have the usual form. If we take, for example, 



the formulation of Lovelock and Rund [|H| , we find that the divergence theorem, a special 
case of Stoke's theorem, takes the following form in an n-dimensional space. If we have a 
vector A = A l ei, the divergence theorem is given in terms of an (n — 1)— form 7Tj defined 
by 

Tij = (-iy+V A • ■ ■ A J' 1 A cu j+1 A • • ■ A u n (2.1) 



Stokes theorem for the form a = -KjA j is 

J da = J a. (2.2) 

G dG 

If we calculate da we find 

da = A)^ 1 A • • • A u n + A^-l^+W A • • • A uj^ 1 A A ■ ■ ■ A w" 

+ • • • + A j (-l) j+ W A • • • A A • • • A du n , 

where A 1 - = e^A*. Using the fact that in the terms containing A> and da/ only the terms 
\C}lj oj k A uji give a non-zero contribution, we find that 

da = (A* ti + A l C l 3l )uj l A • • ■ A u n . (2.3) 

This means that 

J A^u 1 A • • • A uj n = J -A'tfjw 1 A • • • A uj n + J a (2.4) 

G G dG 

It is a moot point whether first term on the right-hand-side of (2.4) should be called a 
"surface term" or not. Here we have chosen to say that only the integral over dG is to be 
regarded as a surface term. 

We would now like to apply the above form of the divergence theorem to Class 
B Bianchi models. If we write the Einstein action for Bianchi models in the basis 
(dr, u 1 , u 2 , u 3 ), where r is a time we choose and {u 1 } is the invariant basis, we find 
that 

1 = ^—1 4 R x P^dr A uo 1 A uj 2 A u 3 , (2.5) 

where 4 g is the determinant of the metric in the above basis. The usual ADM reduction 
gives 

I = ^— ( [tt% - NH± - NiW]dr A lo 1 A lo 2 A u 3 , (2.6) 

167T J 
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where the are treated as an independent momenta, gr^ is the metric on r = const, 
surfaces in the basis {a/},iVj = g 0i , and N = l/(— g 00 ) 1 / 2 . The Hamiltonian and space 
constraints are 

H x = -^g{R + g- l [\{<K k k) 2 - tt%-]}, 

■H' = -2(vrJ + r>^ fe + T%n* k - \%^). (2.8) 

where, as usual, three-dimensional indices are raised and lowered with and gr tJ ' = [g^]" 1 , 
(7 = det[(7jj], comma z means operation with ej, the invariant vector dual to u> 1 , and R is 
the three-dimensional Ricci scalar [T3] , 

<T(r< s ),, - / s (r^), s + /-r*„ri - /t^L - <rc£rt t . (2.9) 

We would like to assume that the Bianchi models are spatially homogeneous, that is that 
all quantities appearing in the action have no space derivatives. Because of the changed 
divergence theorem this is not possible, so we will attempt a "maximal homogenization" 
compatible with Class B models. Spatial derivatives appear in three places: i)n l j , ii) T*^, 
and iii) the spatial derivatives of the that appear in the expresion for T l - k 

1 
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9 % (9lj,k + 9ik,j — 9jk,i) 



+\{C) k - g jP 9 in C p nk - g vk g m Cl 3 ). (2.10) 



We will show that is possible to ignore the in the T* fc and use in the final action 
only T l jk = \{C l j k — gj P g m C v nk — g p kg m Cnj)- From the form of the action (2.7) and of H\ 



it is obvious that the term with 7r l j will cause no problem in any gauge in which iV, = 0, 
and we will always work in such a gauge. All of this means that only the terms in R that 
depend on derivatives of the T l - k will cause problems. It has been shown [IT] that the 



problem with R without derivative terms is that 5(^fg~R) / '5gij ^ —^fg~R % i + \ ^fgg % ^R. We 
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will show that if R is given by (2.9) and we use the divergence theorem (2.5), that indeed 
6(^/gR)/5gij = —^/gR^ + \^[gg % iR. We will use R to denote R written with r* fc in place 
of r i j k . If we set the derivative terms in R equal to zero we find that 

R(T, = 0) = AflrCt^ - l -C a sl C b tJ g st g ab gV 

+\g st 9 pr Cl p C b tT g b3 + Cl r C k tk g rt . (2.11) 

When we vary R we will homogenize after variation, so after variation all derivative terms 
will be put equal to zero. If in (2.9) we replace Y l - k by Y % - k1 we can see that the term 
g rs (r l rl ) iS will not contribute to the variation since T l rl = C l rl and is independent of g^. 
Variation now gives 



Using (2.11) and the definition of T* fci it is not difficult to show that 

5 -^- = lV99' J R-V9R^ (2-13) 

and since in the spatially homogeneous case R = R we see that our assertion that variation 
of R using the proper divergence theorem gives the correct Einstein equations. 

In the following two sections we will apply the above formalism to the simplest classes 
of vacuum models and construct a Hamiltonian formalism. In our case we can follow the 
outline given in Ryan and Shepley where we use the parametrization 



9V = e 2n (e 2 %, (2.14) 
where Q = Q(t) and the matrix is 

)3 = e-^V^V^A^V" 1 ^ 3 (2.15) 



where 





' 


1 


0" 




"0 





0" 




-1 
















1 




_ 





0_ 







-1 


0_ 



(2.16) 



and 

d = diag(/3+ + VSP-,P+ - >/3/3_, -2(3 + ). 
we also define a matrix p^ as 



(2.17) 



Pij = 27r(ef s 7r t s e/) - y (%tt£, 



(2.18) 



where p^ is a function of (p±, p^, p,/,, p#, /3±, 6 1 , 0, ^) such that 

^ — J [p+d(3+ + P-d(3- + p^dip + Ptj>d(f> + pgdO — HdVt]. 



(2.19) 



The Hamitonian H is a function of (p±,p</>,Pip,Pe, P±, 0, 0, ip) equal to 27r(7r^) and is ob- 
tained by solving H± = 0. In those cases where the W are not automatically zero, we add 
the additional constraints TL l = 0. It is not difficult to show that for all Class B models 
TV ^ and 

H 2 = Qp ijPij - 24vr V 6 ^. (2.20) 

The "potential" term now has two parts, R(T, = 0) and g rs (F l rsl ) and the Einstein equa- 
tions are Hamilton's equations in the following form: 



q = 5H/5p, p = -5H/5q, 



(2.21) 



where • = d/dfl and the 5 derivatives include integration by parts by means of the correct 
divergence theorem. An important point is that the term g rs (T l rsl ) can never be converted 
into a pure potential term, because the part of the variation of the form Sg rs (T l rsl ) is set 
equal to zero, while g rs (5T l rs l ) gives a non zero term. There does not exist any pure 



function of which has this property. Note also that a good test of the correctness of 
(2.21) will be (W)' = as an identity. 

In the following two sections we shall study the diagonal (f3 diagonal) and symmetric 
cases (P with one off-diagonal term) to show effect of the term g rs (T l rsl ) on Hamilton's 
equations. 



III. DIAGONAL VACUUM CLASS B MODELS 



It is worthwhile to consider the vacuum diagonal Class B models as a cautionary 
tale, because they are a zoo of pathology. In principle it is easy to write the corrected 
Hamiltonian for the Class B models. We have p^- = p + a\ + p^a 2 , where 





"1 


" 




"V3 





0" 









1 


a 2 = 










(3.1) 







-2 
















The squared Hamiltonian becomes 

H 2 = p\ + p 2 _ + 12ir 2 e- m V([3±) - 2A7r 2 e- 6n g rs (f l rsl ). 



(3.2) 



We use the form of the Ellis-MacCallum [[| scheme for writing the constant Cj fc given in 



Ryan and Shepley [15|. That is 



C) k = e jks m si + 8\aj - 8)a k . 

For all Class B Models a { = (l/2)C l a oc 5f, so only the term g rs (fl s ),3 contributes to H 2 . 
Since m^a,- = by the Jacobi identity, m l3 = for all i. If we define the spin coefficients 



m — m 



of the matrix m , A, B = 1, 2 as a = !(m n + m 22 ), A = 
Hamiltonian (3.2) reduces to 

H 2 = p\ + p 2 _ + 12vr 2 e- 4n e 4/3 +[2(a 2 + A 2 ) cosh (4^3(3. 
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22 \ 



7 = m , the 



+4Aa sinh (4^3(3.) - 2a 2 + 2A 2 + 4 7 2 + 12a 2 ] - 247rV 4 V /3+ [-12a 3 (/3+), 3 - 4v / 37(/3-), 3 ], 

(3.4) 

where we have replaced T l rs by its expression in terms of in the parametrization (2.14) 
and taken the necessary derivatives. Notice that the first of the derivative terms depends 
only on (3 + , and this allows one to integrate by parts. That is, using the correct divergence 
theorem, variation of the part of H containing a 3 gives 

^L(a = A = 7 = 0) = -^{5767r 2 e- 4 °e 4/3 +a 2 - 576vr 2 e - 4 V /3 +a 2 } = 0, (3.5) 

so a reduced Hamiltonian for all diagonal Class B models which gives the equations of 
motion for f}± (i.e. H ^ dH/dQ) is given by 

H 2 = p 2 + + p 2 _ + 127r 2 e~ 4Q e 4/3 + [2(a 2 + A 2 ) cosh(4 v / 3/3_) + 4Aa sinh (4^3/3-)- 

-2a 2 + 2A 2 + 4 7 2 ] + 96v / 3vr 2 e- 4n e 4/3+ 7(/?_) i3 . (3.6) 
This reduction was also noticed by Jantzen [IIJ. No further integration by parts is possible, 



because variation of the remaining derivative term with respect to (3 + must give zero and 
with respect to /3_ must give —192y/3n 2 e~ in e 4:/3+ 'j. 

To study this Hamiltonian we give a, A, 7, and a 3 in Table I. It would seem at this 
point that Types V and VI give simple Hamiltonians because 7 = 0. However, it is 
instructive to apply the test suggested in the previous section, that is, (H 1 )' = 0. For all 
the Class B diagonal models 7i l and H? are identically zero and 7i 3 = reduces to 

6a 3P+ + 2^7^- = (3.7) 
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Table I 

Bianchi Type a A 7 a 3 

III 



1 _ 1 

2 2 



IV i -i -1 

V -1 

vi^_i -1(^ + 1) 



VI Vo -1 



2 2 



It is easy to see that Type V models with a = 7 = are consistent. The reduced 
Hamiltonian is 

H 2 = p 2 + +p 2 _, (3.8) 

and p + = is the space constraint. It is even possible to construct a Hamiltonian that 
obeys H = dH/dCl, that is, p + = implies (3 + = /3+ and a final Hamiltonian is 

H 2 =p 2 _ + 127r 2 e- 4 V /3 +. (3.9) 

This Hamiltonian gives as a general solution for vacuum diagonal Type V models 



p_ = const. 

(3_ = (3°_ + p _n 

H = {p 2 + \2n 2 e- m e A ^} 1 ' 2 = -^-e m ( ^] (3.10) 

127T \dt J 

Let us now examine the Type IV case. If we calculate 6a 3 p + (7 = 0) we find 6a 3 p + = 



6a a 



487r 2 e- 4n e 4/3 +{cosh(4 v / 3/5_) - sinh(4 v / 3/3_)}l ^0 ! 
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This equation is one of the reasons we have qualified the diagonal vacuum Class B 
models as a zoo of pathology. A check of (H 1 )' = shows that models of types III, V, 
and VI^-i are consistent, that is those for which a = A = 0. The two remaining types, 
IV and VII/^o, are examples of another problem in reduced variational principles. We 
have written our variational principle for diagonal metrics by inserting a diagonal metric 
directly into the Einstein action. The problem here is that if the R^, (i ^ j) are not 
identically zero, we will lose the equations — (i ^ j). A rapid calculation shows 
that Rij is diagonal for a diagonal metric for types III, V, and Vlh^-i, while for types 
IV and VII/j^o, R12 7^ 0. This seems to indicate that symmetric metrics (gi 2 7^ 0) will 
always provide consistent vacuum equations if such metrics always have i?i 3 = i? 23 = 
identically. In the following section we show that this is indeed the case, and present a 
consistent Hamiltonian formulation for all symmetric vacuum Class B models. 



IV. THE SYMMETRIC CASE 



In the previous section we have shown that the corrected Hamiltonian as given by 
(2.20) does not always yields equations of motion which are consistent with the space 
constraint equation when one assumes that the metric of the homogeneous hypersurface is 
diagonal. This is certainly the case with diagonal metrics that admit groups of isometries 
of Bianchi types IV and VIL^o- We found out that a common feature appeared in those 
cases where the Hamiltonian equations of motion were not consistent with the space 
constraint. Namely, that some of the non-diagonal components of the Ricci tensor do not 
vanish as one might expect from the diagonality of the space metric. 

In this Section we shall consider symmetric metrics (one nonzero off-diagonal element 

13 



in the metric) and show that all the components of the Ricci tensor for which the corre- 
sponding metric component is zero vanish in this case. Therefore the "pathology" found 
with diagonal metrics is not present in the symmetric case, which tells us that the cor- 
rected Hamiltonian for symmetric metrics must provide the precise equations of motion 
that are consistent with the space constraint equation. This is formally shown by writing 
explicity the corrected Hamiltonian for symmetric metrics and the corresponding Hamil- 
ton equations. Then these equations are used along with the space constraint equation 
to show that the latter is fulflled at all times, and therefore that the equations of motion 
are consistent with the space constraint equation. 

As we are assuming a space symmetric metric, we set g 13 = g 13 = g 23 = 9 23 = and 
then we have only to compute the Ri 3 and R 23 components of the Ricci tensor in the 
order to find out whether they vanish or not. The Ricci tensor is given by 

Rij = -atttj + 2g {j g a k a t + g ik g e sjl m a b 

+g jk g sb a b e sU m lk a b - l/2e sjt e kin m tk m ns 

-^e sjt e kin g sh g ab m tb m na , (4.1) 

with Oj, m 1 - 7 , gij and as defined previously. Then immediately follows that the Ri 3 
component is 

R13 = 93kg 33 m 2k al - ^e 13t e kln m tk m nl - ^e 23t e kln m tk m n2 - 

-^e 13t e kln m tb m na - \^e 23t e kXn g 2k g ab m th m na , (4.2) 

where we have used the fact that a\ = a 2 = (see Ref. [15]), g\ 3 = and the symmetry 
and skewsymmetry of g % i and e^k respectively. Expression (4.2) can be reduced to 

R 13 = l/2e kln g lk g ab m 2b m na , 
14 



and it is straightforward to show that this Ri 3 = 0. 

Now we compute the R 23 component of the Ricci Tensor. By the same arguments 
(viz. ai = a 2 = 0, g 23 = 0, and e^ k = —£ji k ) from (4.1) we obtain 

#23 = 93k9 s3 e s2 im lk al - e 13t e k2n m tk m nl - e 23t e k2n m tk m n2 

-^e 13t e k2n g lk g ab m tb m na - ^e 23t e k2n g 2k g ab m tb m na . (4.4) 
Imposing m 13 = m 31 = m 23 = m 32 = and g 13 = g 23 = 0, (4.4) becomes 

R23 = e k2n m 2k m nl - e k2n m lk m n2 , (4.5) 

or simply R 23 = 0. 

We should point out that the above results, namely i? i3 = and R 23 = 0, hold for all 
symmetric metrics which admit groups of isometries of Class B. 

The form of g^ for the symmetric case is (2.14) and (2.15) with ip — 6 — 0. The 
corrected Hamiltonian for this case is 

H 2 = p 2 + + p 2 + ^L= 24n 2 e- m R, (4.6) 

P+ y ~ sinh(2 v / 3/9-) 

where R, as before, is given by 

R = R + g rs (f l J, h (4.7) 

with the Ricci (curvature) scalar defined by (2.11) and the homogenized connections as 
defined in Sec. II. 

As we did in the diagonal case, one can show that the only contribution from g rs (T l rs ),i 
comes from the term g rs (rf s ) i3 , that is, 
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+23 12 (A2), 3 ]-5 n [5%m 2 1,3 



-9 r2 [9 S % P m 2p is 



(4.8) 



Of course we have used #13 = #23 = g 23 = 0, a\ = a 2 = and the symmetry of m y . 

Now we proceed to obtain the explicit form of (4.8) in terms ot the parametrization 
(2.14 - 2.15). However, in order to have a compact expression for the final result we 
introduce the spin coefficients a, A, 7 of the rotated matrix e^me"^ 3 , where m is given 
by 



m = 



11 12 
m m 



21 22 
m Zi m iA 



We find 



a = hm 11 + m 22 ) 



1 



and 



A = - cos 20(m u — m 22 ) + sin 20m 12 , 



7 = — - sin 20(m n — m 22 ) + cos 20m 12 . 



(4.9) 
(4.10) 

(4.11) 



After a lengthy algebraic computation we end up with the following result for the expres- 
sion (4.8), 

g rs (T l rs ),i = e 2Q e^ + [-12a 3 p + , 3 - 4a/3 7 /?-,3 - 2«0, 3 



+2A sinh(4v / 3/3_)0 i3 + 2a cosh(4 v / 3/3_)0 



3 ■ 



(4.12) 



The scalar curvature R has the same form as in the diagonal case with a, 7, A replaced 
by (4.9 - 4.11), so the final Hamiltonian becomes 



H 2 =p 2 + +p 2 _ + - 



sinh 2 (2 v / 3/9-) 



+ 24vrV 4 V /3+ [(a 2 + A 2 ) cosh(4^3/3_) + 
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+2aA sinh(4V3/3_ 



(a 2 - A 2 ) + 2 7 2 + 6a 2 



-247r 2 e- 4Q e 4/3+ [-12a 3 /3 + ,3 - 4^ 7 /L 3 ■ 
+2A sinh(4V3/3_)0 >3 + 2a cosh(4 v / 3/?_ 



2a0 ; 



(4.13) 



To compute the space constraint equations for this symmetric case we make use of the 



expression (2.8) and impose gis = g = #23 = y = and 7Ti 3 = tx 
Equation (2.8) can be broken down into 



_13 



V"2 3 = 7T 



23 



0. 



- C{^ k = 0, 



(4.14) 



_2fc 



(4.15) 



and 



n 3t 7T k C S — C! 3 TT 3k 

9 ^s^tk ^kj" 







(4.16) 



If one uses the expression in Sec. II for the structure coefficients, then by virtue of the 
symmetric form of g^ and tc^ Equations (4.14) and (4.15) are satisfied automatically. 
However equation (4.16) becomes 



(l 2\ 12 2 12 11 1 1 22 

[Til — 7r|)m — -K\m m + ir^m 



+03(^ + ^-2^) =0. 



(4.17) 



The 7Tj can be obtained, as before, from (2.18). In the Ryan and Shepley [JOd] parametriza- 
tion the trace free part of 71^, i.e. Pij, is given by 



6p 



13 



-4>K A 



[aip + + a 2 p- + a 3 



sinh(2V3/3_) 



with cui and a 2 as before and 



«3 



1 

1 




(4.18) 
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It follows from (4.18) that (4.17) can be written as 

3cf3P+ — Sap^ — 3A coth(2v / 3/9-)p</, + VS'yp- = 0, 

where we have used the spin coefficients previously defined. Now we can check that the 
Hamiltonian (4.13) yields equations of motion which satisfy the space constraint equation 
(4.19) at all times. This we can do by taking the time derivative of (4.19), that is 

3a 3 p+ + 2V3~7p_ - [3a + 3A coth(2v / 3 / 9_)]^ 

+6v / 3Acsch 2 (2v / 3 / 9-)p0/3_ + [67 coth(2V3»4> - 4v / 3>-]0, (4.20) 

where one uses ^ = 27 and ^ = — 2A. Substituting Hamilton's equations as obtained 
from (4.13) for and <j) into (4.20), one verifies that it is indeed zero. 

V. CONCLUSIONS 

In this paper we have investigated the causes why the variational principles for Class 
B Bianchi models breaks down. We have identified as the basic cause of the problem 
that the divergence theorem in a non-coordinated bases does not have the usual form. 
The corrected form of the divergence theorem has been derived, and we have shown 
that applying the usual ADM reduction to the Einstein action for Bianchi models in 
combination with this corrected form of the divergence theorem gives the right Einstein 
field equations for Class B models. 

The above formalism was used to construct a Hamiltonian formalism for diagonal {(5 
diagonal) and symmetric {(5 with one off-diagonal term) vacuum Class B models. 

In the diagonal case we found that the models of Types III, V and VI/^_i are consistent 
in the sense that the space constraint equation is satisfied at all times. However, Types 
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IV and VII/j^o are not. The problem in the latter two is that although a diagonal metric 
is inserted directly into the Einstein action, R 12 is not identically zero and therefore the 
R 12 = equation is lost. 

For the symmetric case we showed that R13 and R23 automatically vanish as a result of 
introducing a symmetric metric. Furthermore, all vacuum Class B models yield equations 
of motion which are consistent with the space constraint equation in the sense previously 
mentioned. 
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